Static dielectric response of icosahedral fullerenes from Ceo to C2160 by an all electron 

density functional theory 
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The static dielectric response of Cqq, Ciso, C240; C540, C720; C96O5 C1500? and C216O fullerenes is 
characterized by an all-electron density-functional method. First, the screened polarizabilities of C60, 
C180, C240, and C540, are determined by the finite-field method using Gaussian basis set containing 
35 basis functions per atom. In the second set of calculations, the unscreened polarizabilities are 
calculated for fullerenes C60 through C216O from the self-consistent Kohn-Sham orbitals and eigen- 
values using the sum-over-states method. The approximate screened polarizabilities, obtained by 
applying a correction determined within linear response theory show excellent agreement with the 
finite-field polarizabilities. The static dipole polarizability per atom in C2160 is ( 4 A 3 ) three times 
larger than that in Ceo ( 1.344 A 3 ). Our results reduce the uncertainty in various theoretical models 
used previously to describe the dielectric response of fullerenes and show that quantum size effects 
in polarizability are significantly smaller than previously thought. 
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I. INTRODUCTION 



Reduction in the spatial dimensions in nanoparticles 
results in a number of interesting properties such as, for 
example, the widening of the band gap, reduction in 
melting temperature, increased magnetic moments etc. 
The fundamental understanding of the size dependence 
of properties of nanoparticles, i.e. at what size the quan- 
tum size effects set in, is important for tailoring these sys- 
tems for possible nanotechnological applications. Stud- 
ies along these lines continue to be frontline research in 
nanoscience 1-9 . 

Fullerenes are hollow cages formed by sp 2 bonded car- 
bon atoms. They are finite analogues of graphene with 
12 pentagonal defects that introduce curvature and make 
them closed structures. Unlike the most nanoparticles 
however, fullerenes are distinctly different in that they 
have all the atoms on the surface whereas solid spher- 
ical nanoparticles have roughly N 1 / 3 atoms on the sur- 
face. This feature results in rapid increase in their vol- 
ume as the number of atoms (i.e. its size) increases. As 
dipole polarizability is proportional to volume, fullerenes 
are most suitable for investigating quantum size effects 
manifested in polarizability. The large dipole polariz- 
ability of carbon fullerenes compared to that of metallic 
or dielectric spheres of same size have been previously 
interpreted as quantum size effect 10 . A number of stud- 
ies on the polarizability of fullerenes have been devoted 
to understanding its size dependence 10-12 . These stud- 
ies have been mainly based on classical models or have 
employed drastic approximations to simplify the polar- 
izability calculations. Unfortunately, the predictions of 
these models are not consistent with each other and show 
rather large variance. For example, the various predic- 




FIG. 1: (Color online) The fullerenes structures optimized at 
ADFT/6-311G** level of theory for which polarizabilities are 
calculated. 



tions of the polarizability of the C2160 (See Table I) are: 
7100 A 3 (Pariser-Parr-Pople model 11 ), 2726 A 3 (bond- 
order-bond polarizability 12 ), 9548 A 3 (Penn model- linear 
response theory 10 ), 17817 A 3 (tight binding 10 ). The pre- 
dictions from different models of fullcrcnc polarizabilities 
and associated quantum size effects differ by almost an 
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order of magnitude. The most reliable methods to de- 
termine electronic structure are the first-principles-based 
quantum-chemical methods. In the preset work, we use 
the first-principles density-functional methods to accu- 
rately determine the first order electric response of car- 
bon fullerenes of sizes from Ceo to C2160 by computing 
their static dipole polarizabilities. These calculations 
provide the most accurate prediction of the dielectric 
response of fullerenes in the size range C60-C2160 and 
show substantial quenching of quantum size effects by 
C 960 in contrast to current understanding. Moreover, the 
fullcrcne polarizabilities computed herein also provide an 
accurate input to the models used in determination of 
the ultraviolet (UV) spectrum for possible detection of 
of hyper-fullcrene particles in space 13 . 



II. THEORY, RESULTS AND DISCUSSION 

The first essential step in obtaining an accurate esti- 
mate of polarizability of fullerene involves accurate de- 
termination of the fullcrcne structure. We optimized the 
geometry of carbon fullerenes using recent fully analytic 
and variational formulation of density functional (ADFT) 
theory 14-19 . Being free of numerical grids, the ADFT im- 
plementation is computationally very efficient and allows 
calculations of matrix elements and energies accurate to 
machine precision. Its functional form is restricted but 
its space of atomic parameters is rich. It has been suc- 
cessfully used to study electronic properties of boron ni- 
trides, aluminum nitride tubes 20-22 , and provides dipole 
moments, total and atomization energies that are com- 
parable to those obtained using some of the most so- 
phisticated density functional models. In this work, we 
optimized the geometries of Cqq, C240, C540, C960, C1500, 
and C2160 fullerenes from the 60 iV 2 family of icosahcdral 
fullerenes. Additionally, two fullerenes (Ciso and C720) 
from the 180A 2 family are also optimized. The struc- 
tures (Cf. 1) are optimized at the all electron level us- 
ing the triple-zeta-quality (6-311G**) orbital basis 23 and 
an exchange-correlation basis set that includes / func- 
tions, and the exchange-correlation parametrization that 
gives an exact Cgo geometry 24 . These approximations 
give accurate geometries in a very efficient manner, which 
we verified by reoptimizing the larger C240 fullerene us- 
ing even larger NRLMOL basis (35 basis functions per 
atom) and using the more complex and popular Perdew- 
Burke-Ernzcrhof (PBE) 25 generalized gradient approxi- 
mation (GGA). Both models give identical C240 geome- 
tries within a root-mean-square value of 10 -3 a.u. for 
the forces. 

For the polarizability calculations the fullerenes struc- 
tures are held rigid i.e., the ionic (vibronic) contribution 
to dipole polarizabilities is ignored. This is a tenable 
assumption as our earlier studies on Ceo and C70 show 
that in fullerenes electronic part determines the total po- 
larizability with negligible contributions from the ionic 
vibrations 26 ' 27 . 



We have determined polarizabilities within two compu- 
tational models. In the first set of calculations polariz- 
abilities of Ceo, Ciso, C240, and C540 are obtained by the 
finite- field method 28-33 , using the NRLMOL code 34-37 . 
In NRLMOL, optimal numerical grids are used to obtain 
the energies and the contributions to matrix elements 
from the exchange-correlation terms. It is well known 
that an accurate description of electric response requires 
use of high quality basis sets. Here, we use large Gaus- 
sian basis sets consisting of 5 s— , 3 p— , and 1 d— type 
Gaussians each contracted from 12 primitive functions 38 . 
This basis set is given in Table II. More details about the 
basis set can be found in Refs.26,38,39. Further inclusion 
of long range d-type basis functions increased polarizabil- 
ity of Ceo by only 0.6%. Therefore, the polarizabilities 
a of Ciso, C240 and C540 were computed without any 
additional diffuse functions to avoid linear dependency 
of basis sets. The total number of basis functions in 
the finite-field polarizability calculations with NRLMOL 
range from 2100 in Ceo to 18900 in C540. The latter is 
comparable to the 38,880 orbital basis functions used in 
our ADFT calculation on C2160- The mean polarizability 
is obtained from the diagonal elements of the polarizabil- 
ity matrix whose components can be determined either 
from the total energies or from the induced dipole mo- 
ments at various field strengths. The polarizabilities ob- 
tained using these two approaches agree within 1%. The 
finite-field a values are given in the third-to-last column 
of Table I. Although the finite-field method is perhaps 
the best approach to calculate polarizabilities, the com- 
putational cost becomes prohibitive for larger fullerenes. 
The cost rises principally due to the very large number of 
self-consistent cycles required to obtain tightly converged 
energies (10~ 9 a.u.) in presence of the external electric 
fields. It is necessary to keep the convergence criteria 
tight to extract meaningful estimates of polarizability 
by the finite-field method. Therefore, we resort to an 
approximate technique that uses the unscreened polar- 
izability obtained by the sum-over-states (SOS) method 
and the linear response theory and random phase approx- 
imation to estimate the screened polarizabilities from the 
&SOS- In the independent electron model, the working 
SOS expression for the xy component of polarizability 
tensor can be written as, 
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Here, the are single particle orbitals and the 

is the single particle energy of the i th orbital. The 
&SOS = ^oixx+cx-yy+cx-zz), is calculated according to Eq. 
(1) once the self-consistent solution of the Kohn-Sham 
equations for each fullerene is obtained using the ADFT 
code. The computational efforts involved in determin- 
ing the large number of transition dipole moments arc 
reduced by making use of group theoretical techniques. 
The calculated otsos f° r fullerenes Ceo through C2160 
are given in the correspondingly labeled column of Ta- 
ble I. The unscreened otsos for Cgo is 292 A 3 , which 
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TABLE I: The average ionic radius r(in A) and polarizabilities (in A 3 ) obtained in various models: bond-order bond polariz- 
ability (BOOP) 12 appp : Pariser-Parr-Pople approach 11 , ap enn : Classical electrostatic within linear response theory (Penn 
model) 10 , TB: tight-binding 10 , asos- Sum-over-states- ADFT (present), ctppA'- Sum-over-states- ADFT-RPA (present), a FF : 
finite- field PBE-GGA (DFT) (present), e: estimates of static dielectric constant for fee structures. 



Fullerene 


r 


BOPP 


appp 




TB 


«sos 


ap.PA 


app e 


C60 


3.55 


76 


80 


64 


81 


292 


79 


82 4.50 


Cl80 


6.13 


227 


209 






1202 


300 


295 4.76 


C24O 


7.07 


303 


306 


343 


581 


1754 


432 


441 4.79 


C54O 


10.27 


681 


866 


1026 


1869 


4823 


1155 


1193 4.71 


C72O 


12.17 


909 


1390 






8000 


1848 


- 4.96 


C96O 


14.03 


1212 


2150 


2185 


4290 


12023 


2745 


- 5.00 


C1500 


17.94 


1893 


4415 






29706 


5690 


- 5.61 


C2I6O 


21.01 


2726 


7100 


9548 


17817 


38288 


8548 


- 5.01 



compares well (282 A 340 and 3 ll 41 ) with previous predic- 
tions with different basis sets and exchange-correlation 
functionals. 

The asos determined using Eq. (1) gives unscreened 
polarizability since the field-induced polarization effects 
(local field effects) are missed in the calculations. Con- 
sequently, the a sos overestimates the experimentally 
accessible screened polarizability 42 . For the smallest 
fullerene (Cgo), asos is roughly 3.6 times larger than 
the screened polarizability in agreement with previous 
reports 40 ' 41 . This overestimation slightly increases (to a 
factor ~ 4) for Ciso, C240, and C540 as can been seen 
from comparison of asos and aFF from Table I. The 
screened polarizability can be deduced from asos using 
the linear response theory 43 and the random phase ap- 
proximation (RPA) by means of the following expression: 



aRPA = 



1 + 



asos 
R 3 



asos- 



(2) 



Here, R is the fullerene radius 43 . For fullerenes (larger 
than Ceo) that are polyhedral in structure, the average 
ionic radius can be used. The average ionic radii r, are 
also given in Table I. Instead of using the ionic radii of 
fullerenes in Eq. (2) which tend to overcorrect the screen- 
ing, it is more appropriate to use an effective radius r+8,. 
The factor 5, accounts for the n electron cloud around the 
bare nuclear skeleton. Its typical values are 1.2-1.3 A 43 . 
Similar methods of determining the screened electronic 
response has been used in past for carbon fullerenes 10 ' 44 , 
multishcll carbon fullerenes- carbon onions 45 , and (cylin- 
drical) carbon nanotubes 43 . However, as we shall see, an 
accurate initial self-consistent treatment is necessary to 
get a meaningful estimate of the screened polarizability 
from Eq. (2). The app a determined using Eq. 2 (with 
S = 1.23A) are listed in Table I. It is apparent that 
aRPA obtained within the ADFT are quite accurate as 
can be gleaned from comparison of uff and app a for 
Ciso? C240 and C540. The app a is within 3% of ajrj? for 
these fullerenes. This excellent agreement provides con- 
fidence in the app a values of the larger fullerenes for 
which a direct calculation of the self-consistent response 
by means of the finite-field method is impractical. 



The electronic structure of these fullerenes show that 
the highest occupied molecular orbital (HOMO) and the 
lowest unoccupied orbital (LUMO) belong to the h u and 
t\ u irreducible representations of the icosahedral point 
group. The lowest HOMO-LUMO (h u —> t lu ) excita- 
tion is forbidden by dipole selection rule and does not 
contribute to the polarizability. The inspection of the re- 
sponse of smaller fullerenes indicate that the low energy 
excitations are strongly screened by application of Eq. 
2. A detailed study of optical spectra will be published 
elsewhere. 

In Table I, the prediction of various models for polar- 
izability are compared. The ap enn is the polarizability 
determined using a semiclassical model that uses classical 
electrostatics and linear response theory. The fullerene 
polarizabilities larger than the ap enn value have been 
interpreted as quantum size effects 10 . The currently ex- 
isting understanding, based on the tight-binding study 10 , 
is that quantum size effects play very important role in 
electronic response of fullerenes up to 03340- In contrast, 
our first principle calculations ( apF and asos in Table 
I) show that by C2160 the semiclassical ap enn is already 
larger than the first-principle value. For the smaller 
fullerenes these effects are also much less pronounced 
than predicted by tight-binding model. For C240, the 
tight binding value of polarizability exceeds finite-field 
DFT value by 32%. Also evident from Table I is that, 
while all models well predict known Cgo polarizability, 
they are inadequate to provide proper description of re- 
sponse of larger fullerenes. The tight binding approxima- 
tion in particular greatly exaggerates the response lead- 
ing to large a values. As same Eq. (2) has been used 
in tight-binding study 10 , it is clear that proper descrip- 
tion of dielectric response require accurate self-consistent 
density functional solution. The table also indicates dra- 
matic increase in polarizability per atom with fullerene 
size. The polarizability per atom systematically increases 
from 1.34 A 3 in C 60 to 4 A 3 in C 2 i6o- There is no ob- 
vious reason why the polarizability per atom will not in- 
crease indefinitely as the rigid structure of fullerenes will 
not allow collapse. Furthermore, the fullerenes become 
energetically more stable with increasing size. The only 
problem is physically making them. 
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The ctRPA values given in Table I are the most accu- 
rate predictions of fullerene polarizabilities. The quan- 
titative agreement with experimental values of polariz- 
abilities require a careful treatment of all possible contri- 
butions to polarizability. The vibrational contributions 
as noted earlier are insignificant 26 . The excellent agree- 
ment between the predicted polarizability and experi- 
mental measurements for Cgo and C70 rule out the role 
of temperature effects 26,27 . Previous studies on carbon 
nanostructures have indicated negative thermal (volume) 
expansion at low temperature 46-48 . We estimate the ef- 
fect of volume expansion by calculating a^pA for all 
fullerenes within 15% of their equilibrium volume. The 
a increases practically linearly with increase in fullerene 
volume within the chosen range of volume expansion (Cf. 
Fig. 3). This, as expected, highlights the importance of 
accurate geometries for accurate prediction of fullerene 
polarizability. 

Finally, using the Clausius-Mossoti relation, we pro- 
vide estimates of the dielectric constant of crystalline 
solids of these fullerenes. Assuming that like C@o, all 
larger fullerenes form an FCC crystal, the number den- 
sity of fullerene is then 4/V, where V=(-\/(8)r) 3 . The 
fullerene radius used here is the effective radius as used 
above (to include the 7r-electron spill-out effect). Using 
a value of S = 1.23 A results in a smaller lattice con- 
stant than the experimental value (14.17 A) 49 but gives 
a dielectric constant of 4.50, which is in good agreement 
with experimental values (4.0 - 4.5) 50 . The dielectric con- 
stants calculated for the larger fullerenes using the same 
S value are in the range 4.7-5.6 (See last column in Table 
I). These values may vary slightly with the 6 but the 
trend should be the same. 



III. CONCLUSIONS 

In summary, static electric response of the icosahe- 
dral Ceo, C180, C240, C540, C720, Cg6o, C1500, and C216O 
fullerenes is studied in detail by an all-electron first 
principles density-functional methods. Quantitative es- 
timates of the dipole polarizability determined herein 
provide better understanding of electronic response of 
fullerenes amongst earlier inconsistent predictions from 
various models. This work shows that quantum size ef- 
fects in polarizability are substantially quenched by C2160 
in sharp contrast to previous tight-binding predictions. 
These calculations also signals the beginning of an era 
where tens of thousands of basis functions will be used 
for accurate electronic structure calculations. 
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volume, r 3 , as a function of fullerene radius r. The squares 
represent polarizability obtained by the finite-field method. 



computer time at the UTEP Cray acquired using ONR 
05PR07548-00 grant. 



1 H. Haberland, T. Hippler, J. Donges, O. Kostko, 
M. Schmidt, and B. von Issendorff, Physical Review Let- 
ters 94, 035701 (pages 4) (2005). 

2 L. E. Ramos, J. Furthmiiller, and F. Bechstedt, Physical 
Review B (Condensed Matter and Materials Physics) 72, 



045351 (pages 8) (2005). 
3 X. Wang, R. Q. Zhang, S. T. Lee, T. A. Niehaus, 
and T. Frauenheim, Applied Physics Letters 90, 123116 
(pages 3) (2007), URL http : //link . aip . org/link/?APL/ 
90/123116/1. 



5 



to 

N 
's— 

JS 
o 

Q. 



to 

0} 




5 10 

% increase in volume 

FIG. 3: (Color online) The percent increase in polarizability, 
otRPA, as a function of percent increase in volume, r , for 
Cgeo fullerene. 



4 L. Yang, C. D. Spataru, S. G. Louie, and M. Y. Chou, 
Physical Review B (Condensed Matter and Materials 
Physics) 75, 201304 (pages 4) (2007). 

5 J.-Y. Raty, G. Galli, C. Bostedt, T. W. van Buuren, 
and L. J. Terminello, Physical Review Letters 90, 037401 
(pages 4) (2003), URL http://link.aps.org/abstract/ 
PRL/v90/e037401. 

6 N. D. Drummond, A. J. Williamson, R. J. Needs, and 
G. Galli, Physical Review Letters 95, 096801 (pages 4) 
(2005), URL http://link.aps.org/abstract/PRL/v95/ 
e096801. 

7 Y. Volokitin, J. Sinzig, L. J. de Jongh, G. Schmid, M. N. 
Vargaftik, and I. I. Moiseevi, Nature (London) 384, 621 
(1996). 

8 G. A. Breaux, C. M. Neal, B. Cao, and M. F. Jarrold, 
Physical Review Letters 94, 173401 (pages 4) (2005), URL 
http: //link. aps . org/abstract/PRL/v94/el73401. 

9 I. Vasiliev, S. Ogut, and J. R. Chelikowsky, Physical Re- 
view Letters 86, 1813 (2001), URL http: //link. aps. org/ 
abstract/PRL/v86/pl813. 

10 G. K. Gueorguiev, J. M. Pacheco, and D. Tomanek, Phys. 
Rev. Lett. 92, 215501 (pages 4) (2004). 

11 A. Ruiz, J. Breton, and J. M. G. Llorente, Journal of 
Chemical Physics 114, 1272 (2001). 

12 Y. H. Hu and E. Ruckenstein, The Journal of Chemical 
Physics 123, 214708 (pages 4) (2005). 

13 S. I.Groth, The Astrophysical Journal 661, L167 (2007). 

14 K. S. Werpetinski and M. Cook, Phys. Rev. A 52, R3397 
(1995). 

15 K. S. Werpetinski and M. Cook, The Journal of Chemi- 
cal Physics 106, 7124 (1997), URL http: //link, aip.org/ 
link/?JCP/106/7124/l. 

16 B. I. Dunlap, Journal of Physical Chemistry a 107, 10082 
(2003), pT: J. 

17 R. R. Zope and B. I. Dunlap, Physical Review B 71, 193104 
(2005), pT: J. 

18 R. R. Zope and B. I. Dunlap, Journal of Chemical Theory 
and Computation 1, 1193 (2005), pT: J. 

19 R. R. Zope and B. I. Dunlap, Journal of Chemical Physics 
124, 044107 (2006), pT: J. 

20 R. R. Zope and B. I. Dunlap, Physical Review B 72, 45439 



(2005) . 

21 R. R. Zope, T. Baruah, M. R. Pederson, and B. I. Dunlap, 
Arxiv preprint physics/0407031 (2004). 

22 R. R. Zope and B. I. Dunlap, Chemical Physics Letters 
386, 403 (2004), pT: J. 

23 A. D. McLean and G. S. Chandler, The Journal of Chemi- 
cal Physics 72, 5639 (1980), URL http://link.aip.org/ 
link/?JCP/72/5639/l. 

24 B. I. Dunlap and R. R. Zope, Chemical Physics Letters 
422, 451 (2006), pT: J. 

25 J. P. Perdew, K. Burke, and M. Ernzerhof, Physical Review 
Letters 77, 3865 (1996), pT: J. 

26 M. Pederson, T. Baruah, P. Allen, and C. Schmidt, Journal 
of Chemical Theory and Computation 1, 590 (2005), ISSN 
1549-9618. 

27 R. R. Zope, Journal of Physics B: Atomic, Molecular and 
Optical Physics 40, 3491 (2007), URL http : //stacks . 
iop . org/0953-4075/40/3491. 

28 H. D. Cohen and C. C. J. Roothaan, The Journal of Chem- 
ical Physics 43, S34 (1965). 

29 R. J. Bartlett and G. D. Purvis, Phys. Rev. A 20, 1313 
(1979). 

30 H. A. Kurtz, J. J. P. Stewart, and K. M. Dieter, J. Comput. 
Chem. 11 (1982). 

31 J. Guan, M. E. Casida, A. M. Koster, and D. R. Salahub, 
Phys. Rev. B 52, 2184 (1995). 

32 S. A. Blundell, C. Guet, and R. R. Zope, Phys. Rev. Lett. 
84, 4826 (2000). 

33 K. Jackson, M. Pederson, C.-Z. Wang, and K.-M. Ho, 
Phys. Rev. A 59, 3685 (1999). 

34 M. R. Pederson and K. A. Jackson, PHYSICAL REVIEW 
B. CONDENSED MATTER 43, 7312 (1991), pUBM: 
Print; ppublish. 

35 M. R. Pederson and K. A. Jackson, PHYSICAL REVIEW 
B. CONDENSED MATTER 41, 7453 (1990), pUBM: 
Print; ppublish. 

36 K. Jackson and M. R. Pederson, PHYSICAL REVIEW 
B. CONDENSED MATTER 42, 3276 (1990), pUBM: 
Print; ppublish. 

37 K. Jackson and M. R. Pederson, Physical Review B 42, 
3276 (1990), pT: J. 

38 D. Porezag and M. R. Pederson, Physical Review a 60, 
2840 (1999), pT: J. 

39 M. R. Pederson and T. Baruah, Lecture Series in Computer 
and Computational Sciences 3, 156 (2005). 

40 E. Westin, A. Rosen, G. T. Velde, and E. J. Baerends, 
Journal of Physics B: Atomic, Molecular and Optical 
Physics 29, 5087 (1996). 

41 M. R. Pederson and A. A. Quong, Phys. Rev. B 46, 13584 
(1992). 

42 R. R. Zope, T. Baruah, M. R. Pederson, and B. I. Dunlap 
(2007), cond-mat/0701466. 

43 L. X. Benedict, S. G. Louie, and M. L. Cohen, Phys. Rev. 
B 52, 8541 (1995). 

44 G. F. Bertsch, A. Bulgac, D. Tomanek, and Y. Wang, 
Phys. Rev. Lett. 67, 2690 (1991). 

45 A. Ruiz, J. Breton, and J. M. G. Llorente, Physical Review 
Letters 94, 105501 (pages 4) (2005). 

46 S. Brown, J. Cao, J. L. Musfeldt, N. Dragoe, F. Cimpoesu, 
S. Ito, H. Takagi, and R. J. Cross, Physical Review B (Con- 
densed Matter and Materials Physics) 73, 125446 (pages 6) 

(2006) . 

47 P. K. Schelling and P. Keblinski, Phys. Rev. B 68, 035425 
(2003). 



6 



M. Yoshida, Y. Kurui, Y. Oshima, and K. Takayanagi, Jap. 
J. Appl. Phys. 46, L67 (2007). 

P. C. Eklund, A. M. Rao, and Y. W. et al., Thin Solid 
Films 257, 7590 (1995). 



A. F. Hebard, R. C. Haddon, R. M. Fleming, and A. R. 
Kortan, Applied Physics Letters 59, 2109 (1991). 



